基于最大球算法定量表征土体空间孔隙网络     

解朝阳  张巍  姚东方  王坚  朱雨辰 《工程地质学报》2020,28(1):60-68

孔隙网络控制着土体渗流、排水固结与基质吸力等重要工程性质。本文介绍了多孔介质孔隙网络最大球建模基本原理与算法。以显微CT扫描振捣干法生成的南京粉砂试样为例,采用最大球算法建立了试样三维重构模型表征单元体(REV)的空间孔隙网络球棍模型,计算得到了样品REV尺度的孔隙网络参数,统计发现,孔隙半径、喉道半径、孔隙配位数、孔隙截面形状因子、喉道截面形状因子与喉道长度等孔隙参数均近似服从正态分布,孔隙体积近似服从衰减型指数分布。孔隙与喉道半径分别分布在100μm与65μm以内,两者数学期望分别为40. 0μm与18. 0μm;配位数分布在25以内,数学期望为5. 1;孔隙与喉道截面形状因子分别分布在0. 01~0. 04与0. 01~0. 05的区间内,两者数学期望分别为0. 019与0. 033;喉道长度分布在100~800μm以内,数学期望为292. 22μm。同时发现,样品中体积小于1. 5×10^7μm^3的小孔隙数量超过90%。本方法可应用于土体细观孔隙结构的定量表征。  相似文献   

Crack nucleation in saturated porous media     

Luc Dormieux  Eric Lemarchand  Thomas Carlioz  Laurent Jeannin 《国际地质力学数值与分析法杂志》2020,44(5):553-574

Crack nucleation has been the subject of important contributions in the last two last decades. However, it seems that few attention has been granted to the case of saturated porous media. This is the question addressed in the present paper which is devoted to nucleation in traction mode. From a physical point of view, nucleation is a sudden phenomenon, so that the material response is both adiabatic and undrained. In the spirit of the variational approach, the nucleated crack is viewed as the final state of a region of space in which the material undergoes a full damage process. In traction mode, the opening of a saturated crack in undrained condition induces a drop of fluid pressure. In case of low fluid compressibility, the presence of the fluid delays the brittle failure usually associated with nucleation, as long as the fluid pressure remains above the saturation vapor pressure. Nucleation is therefore possible only if a partial vaporization of the fluid takes place.  相似文献   

Clogging of saturated porous media by silt‐sized suspended solids under varying physical conditions during managed aquifer recharge          下载免费PDF全文

Xinqiang Du  Xueyan Ye  Xiaowan Zhang 《水文研究》2018,32(14):2254-2262

Managed aquifer recharge is an effective method for utilizing excess flood flows, but clogging of porous media is a limiting factor in the implementation of this water storage technique. In recent years, much research on the physical clogging of porous media during artificial recharge has been conducted. However, the understanding of clogging due to silt‐sized suspended solids (SS) is still inadequate, especially under varying physical conditions. Here, we subjected sand columns to controlled rates of flow and SS suspensions to investigate the influence of media size, SS size, SS concentration, and flow velocity on the clogging of porous media by silt‐sized SS. The results show that the diameter ratio of SS particles to sand grains is the dominant factor influencing the position of physical clogging. As pore velocity increased, the mobility of silt‐sized SS was enhanced and retention in the porous media decreased noticeably. The spatial retention profiles in the porous media were found to vary greatly at different flow velocities. The SS concentration of the infiltrating suspension also dramatically influenced the mobility and deposition of silt‐sized SS particles, such that high concentrations accelerated the clogging process. As the different physical factors changed, the breakthrough curves and retention profiles of silt‐sized SS particles changed obviously and the mechanisms of retention differed. On the whole, clogging position is mainly determined by particle size ratio, but clogging rate is dominated by a variety of factors including particle size ratio, SS concentration, and flow velocity.  相似文献   

含孔隙介质的分层半空间表面瑞利波的衰减特性          下载免费PDF全文

阎守国  谢馥励  张碧星 《地球物理学报》2018,61(2):781-791

以分层半空间内部含有一层孔隙介质为物理模型进行数值计算，研究半空间表面瑞利波的传播和衰减特性.为更加接近实际，结合瑞利波的激发特性，确定了瑞利波的主衰减曲线，并主要以此进行规律分析.针对速度递增和含低速层这两种典型的地质模型，讨论了瑞利波的传播衰减在不同地质模型下的特性，并分析了各自的规律.结果表明，在这两种模型下瑞利波的主衰减曲线都受孔隙介质所处空间位置影响产生比较明显的变化，但衰减系数极大值对应的波长与模型的表层厚度存在较明显的线性对应关系，利用这一关系，可以在实际勘探中快速得到表层介质厚度.另外，通过对比分析还可以看到，瑞利波主衰减曲线随孔隙介质的孔隙度和渗透率的变化都强于主频散曲线的变化，表明衰减曲线对孔隙度和渗透率的变化更加敏感，理论上更加适合进行介质参数反演工作.综合对比结果，我们认为瑞利波主衰减曲线中包含了更丰富的介质参数信息，如果能够有效利用，将可以提高瑞利波勘探的准确性和应用范围.  相似文献   

含流体致密砂岩的纵波频散及衰减:基于双重双重孔隙结构模型描述的特征分析          下载免费PDF全文

郭梦秋  巴晶  马汝鹏  陈天胜  张琳  庞孟强  谢剑勇 《地球物理学报》2018,61(3):1053-1068

致密砂岩普遍具有低孔、低渗及微裂缝发育的地质特征，并且呈现出很强的非均匀性.致密砂岩储层与常规砂岩储层比较，具有明显的岩石物理性质、渗流力学性质方面的差异.致密砂岩内部的非均匀性对弹性波频散、衰减有显著影响，其中包括孔隙结构的非均匀性，即岩石内部孔隙参数的不均一性，以及孔隙内部不相混溶流体的非均匀分布；此外，非均匀性的尺度也决定了波出现显著频散与衰减的频段.综合考虑致密砂岩孔隙结构非均匀性及流体斑块状饱和的非均匀性，本文采用双双重孔隙介质结构模拟了致密砂岩的弹性波响应，分析了同时具备两类非均质性岩石中的波传播特征.调查分析了两组分别来自中国鄂尔多斯盆地苏里格气田及四川盆地广安气田的不同类型致密砂岩储层的岩芯超声波实验数据，给出了岩石样本的弹性波速度频散与衰减曲线.结果显示理论模型预测结果与完全饱和、部分饱和岩石的实验数据吻合良好.对两个地区致密砂岩岩芯数据进行对比分析，苏里格致密砂岩样本总体上比广安致密砂岩渗透率高，在各孔隙度范围内，特征模拟显示苏里格样本的裂隙尺寸明显大于广安样本.广安致密砂岩在低孔隙度范围内发育了更多、更小的颗粒裂隙/接触.致密砂岩的速度频散与衰减结果受流体黏度、晶体破裂及流体斑块状饱和的共同影响.此外，孔隙度越大，部分饱和岩石中斑块状饱和机制对总衰减的贡献越低，与之相对，结构非均质性所占的比重则有所增强.  相似文献   

非线性波消波及波场分布研究   总被引：1，自引：1，他引：0  

王巧莎  李明海 《海洋科学》2018,42(6):30-38

本文基于雷诺平均N-S方程,并结合RNG k-e方程建立了粘性数值波浪水槽,对不同波陡、不同相对水深、不同相对波高的非线性规则波的阻尼消波问题和波场分布进行研究。文中提出了两种描述消波区内部阻尼变化的阻尼函数,分别适用于小波陡情形和高波陡情形。研究结果表明,小波陡组消波区可设为一个波长,阻尼系数取10~4～10~5即可满足消波要求,计算结果与实验结果及造波理论吻合良好;高波陡组消波区可设为两个波长,阻尼系数取10~4～10~5亦可满足消波要求,计算结果与实验结果吻合良好。此外,当波陡较小时,波场内反射情况的小幅改变即可对整个波场造成影响,特别是当水深较浅时这种影响极为明显,需谨慎考虑。当波陡较大时,水波能量较高,整个波场沿水波传播方向可观测到明显的衰减现象,在具体试验中需进行考虑。  相似文献   

基于纵横波保幅分离的粘弹介质弹性波正演模拟     

侯志强  尹文笋  李键  孙永壮  刘云 《中国海洋大学学报(自然科学版)》2020,(1):82-92

地下介质往往表现为粘弹特性,研究基于粘弹假设的数值模拟方法对于正确认识地震波的传播规律和提高地震勘探精度具有重要意义。波动方程正演过程中纵横波的保幅解耦是研究准确的粘弹介质中地震波传播机理的前提,基于散度和旋度算子的纵横波解耦方法会使波场的相位和振幅产生畸变,且解耦后的波场在极性反转位置上无法与分离前混合波场各分量对应。在散度和旋度算子上再做一次梯度和旋度运算的波场分离方法虽然能够克服上述缺陷,但存在保幅性差等问题。本文从粘弹介质中的一阶速度-应力方程出发,推导了矢量纵横波分离的波数域表达式,结合有限差分思路给出了其在空间域的求解方法。本文方法利用纵横波的传播速度对现有的矢量波场分离方法进行振幅校正,并将校正结果分别作为纵波与横波对时间的二阶偏导,实现了粘弹介质中的纵横波分离。模型试算结果表明,本文方法能够克服现有方法的缺陷,获得更具保幅性的波场分离结果。  相似文献   

饱和土体-地下综合管廊结构地震响应分析          下载免费PDF全文

蒋录珍  李双飞  石文倩  安军海  陈艳华 《地震工程学报》2018,40(6):1224-1230

将土体视为固-液两相介质,基于饱和土体有效应力原理,建立饱和土体-地下综合管廊结构体系相互作用动力模型:在地应力平衡的静力状态下采用Duncan-Chang非线性弹性本构模型,在地震波作用的动力状态下采用Davidenkov非线性黏弹性本构模型;考虑饱和土体黏弹性动力人工边界条件,将地震动作用转化为作用在人工边界节点上的动力荷载。模型考察不同地震波时程、地震波加速度峰值、入射角度、孔隙率以及地应力场的影响,得出如下结论:(1)地震波的卓越周期与场地卓越周期相近时引起结构上的变形最大;随着地震波加速度峰值的增大结构变形增大;随着地震波入射角度的增加结构变形增大,地震波斜入射情况下产生的行波效应使得结构变形最大。(2)土体材料的孔隙水压力是影响地震中结构变形的主要因素之一。(3)将土体材料考虑为单相介质时结构上的变形要比考虑为固-液两相介质时大得多,直接将饱和土体场地中得到的地震波等效荷载施加到单相土介质-结构动力相互作用模型上,能够得到与完全基于有效应力法一致的结果。  相似文献   

一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

A semi-empirical approach to model pressure dependence of elastic moduli in granular media accounting for variations of coordination-number and Poisson-ratio     

Angelo Sajeva  Davide Scarpellini  Pamela Tempone  Simone Capaccioli 《Geophysical Prospecting》2019,67(4):872-887

The effective medium theory based on the Hertz–Mindlin contact law is the most popular theory to relate dynamic elastic moduli (or elastic velocities) and confining pressure in dry granular media. However, many experimental results proved that the effective medium theory predicts pressure trends lower than experimental ones and over-predicts the shear modulus. To mitigate these mispredictions, several evolutions of the effective medium theory have been presented in the literature. Among these, the model named modified grain contact theory is an empirical approach in which three parametric curves are included in the effective medium theory model. Fitting the parameters of these curves permits to adjust the pressure trends of the Poisson ratio and the bulk modulus. In this paper, we present two variations of the modified grain contact theory model. First, we propose a minor modification in the fitting function for the porosity dependence of the calibration parameters that accounts for non-linearity in the vicinity of the critical porosity. Second, we propose a major modification that reduces the three-step modified grain contact theory model to a two-step model, by skipping the calibration parameter–porosity fit in the model and directly modelling the calibration parameter–pressure relation. In addition to an increased simplicity (the fitting parameters are reduced from 10 to 6), avoiding the porosity fit permits us to apply the model to laboratory data that are not provided with accurate porosity measurements. For this second model, we also estimate the uncertainty of the fitting parameters and the elastic velocities. We tested this model on dry core measurements from literature and we verified that it returns elastic velocity trends as good as the original modified grain contact theory model with a reduced number of fitting parameters. Possible developments of the new model to add predictive power are also discussed.  相似文献